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Abstract 


This report gives the results of a number of attempts to find under- 
lying principles uniting the many methods for generating Liapunov functions 
for nonlinear systems of ordinary differential equations. Our conclusions 
indicate that if general principles exist for the generation of Liapunov 
functions they lie very deep s and are not to be found among the methods 
now in use. , 


iii 


a 



Summary 


The report falls naturally into two parts. Part I gives a critique 
of the methods now in use for the generation of Liapunov functions for 
nonlinear systems of ordinary differential equations. The critique is 
illustrated by a technique for generating Liapunov functions which we call 
the "Conservative Spring Concept". This technique, which we emphasize 
is not really a new method, is capable of solving a large number of 
previously published examples. 

Part II of this report is a description of some partial results 
obtained for the stability analysis of systems of autonomous ordinary 
differential equations whose right hand sides are homogeneous polynomial 
forms . 

The report is concerned exclusively with the problem of Liapunov 
stability or asymptotic stability of an isolated singular point, usually 
the origin. 
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I, The purpose of this research was to try to uncover any general principles 
that may exist for the systematic generation of Liapunov functions. After 
looking for quite some time at the various methods used to generate Liapunov 
functions for the analysis of the stability of systems of autonomous differ- 
ential equations it is our opinion that, despite the great bulk of published 
material, little progress has been made. Our efforts to deal with the 
problem of the systematic generation of Liapunov functions has led us to feel 
that a specific, rather than a general, approach is most likely to succeed 
in the long run. It is, in fact, the very generality, or seeming generality, 
of the various methods now in use for finding Liapunov functions which is 
their principal weakness. The authors of the various methods claim general- 
ity, but they always end up applying their methods to basically the same 
classes of second, third, and a few fourth order systems.* The exceptions 
to this rule, which we know of, are high order systems describing certain 
rigid body motions, and control theory problems with high order linear parts 
combined with a single nonlinear element. The methods of obtaining Liapunov 
functions in these two areas are for the former, to consider a linear combina- 
tion of first integrals, and for the latter to consider a quadratic form plus 
an integral for the nonlinear element. 

Most of the methods studied are, as claimed, very general in .principle, 
but they all suffer from the following defects in practice. 

(i) A large amount of computation is required for high order 
systems. This is the least serious defect. 

*It should be understood that we are here talking of practical procedures 
and not of the very substantial body of theoretical results. 
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(ii) At various stages of the application it is necessary to 
check on the sign definiteness of one or more functions. For functions 
other than quadratic forms this is an unsolved problem presenting 
formidable difficulties. 

(iii) At some stage of every method studied the practitioner 
must make an apt, often ingenious, choice of one or more functions. 

This is the chief difficulty, which we elaborate on in the next 
paragraph. 

Most methods for the generation of Liapunov functions are 
motivated by some analogy with concepts which work well for simple and/or 
special types of systems. For example, the variable gradient methods are 
motivated by the use of quadratic forms which work well for linear systems; 
methods which utilize a first integral (either of the system studied or a 
modification of that system) are motivated by a generalization of the total 
energy of the system. In each case the analogy works for a while, but 
leads to ever increasing complexities as the problem to which it is applied 
gets further removed from the type of system which gave rise to the technique. 

We illustrate the above idea by presenting a."new method" of generating 
Liapunov functions which we call the Conservative Spring Concept of Generating 
Liapunov Functions because of its initial motivation. The method, developed 
by Dr. R. Drake, appears on the surface to have great generality s in ce by 
its use a large number of interesting examples presented by the authors of a 
variety of methods are worked out, usually with far less difficulty than was 
presented by the original method. However, we believe that this is due 
mainly to the fact that the examples are either of low order or contain only 
a single simple nonlinearity. Our examples are chosen so as to illustrate 
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the evolution of the method. Each new type of equation is further removed 
from the original idea. As the analogy grows weaker additional assumptions 
must be introduced until finally there are gaps to be filled in by a judi- 
cious choice of functions. 

In summary, we believe at this point that effective progress in the 
problem of generating Liapunov functions has not been made mainly because 
the proposed methods are too general. We feel .that one should consider 
each order of nonlinear system in turn with respect to classifying the 
various types of nonlinearities according to how they affect the stability 
of the system. The Conservative Spring method presented in this report is 
not really a new method. It serves to illustrate the pitfalls of general- 
ity and also serves to illustrate the possibilities of the above type of 


classification. 



• - The Conservative Spring Concept of Generating Liapunov Functions 

1. Motivation. 

The differential equation governing the motion of an undamped nonlinear 
spring is • 

, j " i 

(1) x + g(x) = 0 j g(0) = 0 . • 

The total energy of the system at any time is given by the function 

* 2 ; 

' V = | + G(x) , . ‘ 

• x • • , 

where G(x) = J" g(s) ds. Differentiating V with respect to time along 

o ■ 

a trajectory of (1) yields ' 

• * •• * • • 

V. = x x + g(x)x = x (-g(x)) + X g(x) = 0 . 

Hence if we require that 

G(x) > 0 for x =f= 0 , 

the function V becomes a Liapunov function for (1) and we may conclude that 

(1) is stable. 

2. Lienard’s equation. 

We may consider Lie'nard’s equation 

•• • 4 

(2) x + f(x) x + g(x) = 0 

as representing the motion of a non-conservative spring. We attempt to cast 


Let 


equation (2) into a form which resembles (1), 


-f 


F(x) = J f (s) ds 

o 


-i 


G(x) = J g(s) ds , 
o 


and rewrite (2) as 


d_ 

dt 




+ g(x) 


0 


By analogy with (1) we choose the function 


(3) 


V 




+ G(x) 



as a candidate for a Liapunov function. Note that now V no longer represents 
the total energy. Differentiating V with respect to t along a trajectory of 
(2) yields 


V = -g(x) F (x) . 

We may conclude that (2) is asymptotically stable in a region defined by 
the conditions 

g(x) F(x) > 0 for j x | < a, x ={= 0 

G(x) < H | x | < a 

S2 = { (x, x) / V < £ } 

Complete details of the analysis may be found in f 5 ] . 
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We now compare, by means of several examples, the proposed Liapunov 
function given in equation (3) with the Liapunov function obtained by other 
methods.* 

* . c • 

Example A. Van-der Pol's equation 


2 * 

x + e (1-x ) x +, x = 0 

x 3 

F(x) = £ ^ (1-s 2 ) ds = e (x - ~ ) 

o . 

r x 2 

G(x) = J sds = 2 * 

o 

' V* 3-122 

V = f [x + e " f )J + f . . . 

• 2 x^ 

V = - g(x) F (x) = - ex (1 - -j ) . 

The Liapunov function given by this analysis is equal to 1/2 of the 
Liapunov function given by Infante ^2;p 74}. 

Example B. Infante {2;p 78} considers a non-symmetrical oscillator 
with equation 

• 2 ' 
x + ax + bx + X = 0 ; a, b constants > 0. 


Here 

■nr \ nr \ -b 2 , 1 3 

F(x) = ax, G(x) = — x + x 

In accordance with (3)‘ we take 


* The method with which we are comparing will be denoted by the name of its 
author; ' references are to section and page of the report [6 3 . 
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v “ I [* + 8X ] 


2 J b 2 1 3 

+ 2 * + 3 x 


' 2 2 
V = - (bx + x ) ax = - ax (b + x) 


This is the same Liapunov function as obtained by Infante. 

Example C. Infante (2;p 81 } considers a nonlinear damped pendulum 
with equation 

•• • 

x + (e cos x) x + sin x = 0 ; e > 0 . 

Here 

F(x) = £ sin x , G(x) = 1 - cos x 

and 


V 


= ~ J^x + e sin xj 


+ (1 - cos x) , 


• 2 

V = - e sin x . 


This is the same Liapunov function as obtained by Infante. 
Example D. A globally stable oscillator. 


Here 


x + e (1 - x^ .+ x^) x + x^ — 0 


3 5 4 

F(x) = £ (x - ~ ) , G(x) = — 


and 


* |[x + e (X - 


3 5-2 

5- + *_ > 

3 5 


* 4 ,, x . 5 v 

V = -£X (1 - -r~ + x ) 

5 

Again, our Liapunov function is the same as that given by Infante (2;p 82} 
Szego {7;p 46} using a method not related to the energy concept obtains th 
function 
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rj 1 4 ’ 2 

v s = - X + X , 


' '2 2 4 

V = 2 ex (1 - x + x ) . 
s 


Example E. Schultz and Gibson {5;p 28} consider the following example 
by means of the variable gradient technique. 


x + (1 + f(x) + x f(x)) x + 8 x f(x) = 0 


Here 


F(x) = x + y f(s) ds + ^ s f(s) ds = x (1 + f(x)) 
o o 


i 


X 


G(x) = 8 J s f(s) ds . 

o 


In accordance with (3) we set 

V = [x + x(l + f (x))J 2 + 8 J s f(s) 


x 


ds 


V = - 8 x f (x) (1 + f (x) ) . 


The conditions for global asymptotic stability are 
x f (x) >0 ; x =|= 0 

8 > 0 

|f (x) | >_ 1 

The Liapunov function obtained by Schultz and Gibson is 



\2 


W = (x + x) + 2 C s [f(s) + s ft s )] ds + 2g j"' s f(s) ds 

J o o 


W 


= - 2 £(f(x) + x f'(x)) x 2 + 3 x 2 f(x)J 


The Liapunov function W gives the following conditions for global asymptotic 
stability 

• g > 0 ■ ' ... 

f (x) + xftx) _> 0 
f (x) > 0 for x =[= 0 

3. A more general second order equation. 

The Equation 

(4) x + f (x) x 2 + g(x) x + h(x) = 0 

is one step removed from the Lienard equation; we therefore aim to put it 

into a form similar to (2). 

_ x 

Let y(x) = exp J f(s) ds J = exp [j(x)J 

o . 

Multiplying (4) by y gives 

(5) (x + f(x) x 2 ) y(x) + g(x) y (x) x + h(x) y(x) = 0 


Define 


x x 

G (x) = J g (s) y (s) ds , H (x) = J h(s) y(s) ds , 
o o 
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Since y = x f(x) y(x) we may write (5) in the form 


dt 


[iy + S] 


x y + GJ + h(x) y(x) = 0 


By analogy with (3) we take as a candidate for a Liapunov function 


v = | [_ * y ^ x - > + G ^J 2 + 

Differentiating with respect to t along a trajectory of (4) yields 
V = x h(x) ,^exp [F(x)] - exp [2F(x)]^ - h(x) G(x) exp [F(x)3 . 


V would have the same form as in (3) if the first term above were zero. 
This can be accomplished by modifying V by use of the function 


H(x) 


= J" h(s) exp [2F(s)3 ds 
o 


in place of H(x) . 

Our Liapunov function now becomes 
(6) V = — £ x exp [F(x) ] + G(x)j + H(x) 


V = - h(x) G(x) exp [F(x)] 

‘ A 

The function G(x) is not really needed in (6); if we delete it from the 
definition of V in (6) we get the alternate form 
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(7) 


V = x 2 exp [2F (x) ] + H(x) s 


- V = - x 2 g(x) exp [2F(x) ] . 

Example, Infante {2;p 83} considered the problem of liquid motion in 
a surge tank, and derived the governing differential equation 


x + 


or 

8 


*2 . 8 T 2a 2 . , . 2a 2 

x + Tfi+i* > L"T - i + — ^ i x 


8 


0* 


c«- 


8 + 


8 ' 


il- 


(i+x): 

Using the alternate form (7) we obtain the function 
V = \ x 2 exp [ 2 x] + J [ s - 8 ^-^^exp [ 2 -J- sj ds, 

O _n ' 


V = - 


‘2 

8 x 


f 2a 2 


2 

i 0 a * 

1 

L 8 

- 1+ T X J 

exp [2 — 

"J 


for the analysis of the stability of (8) , which is the same as the 

Liapunov function given by Infante. 

4. The general second order equation. 

We ■ turn now to the most general second order equation, 

•• • * * 

(8). x.+ X(x,x) = 0 ; X(0,0) = 0 . 


- 12 - 



In order to get equation (9) into a form amenable to our. previous analysis 
9X 9^X 

we assume that and „ exist and are continuous. Write. 


ax 


3x 2 


X(x,x) = h(x) + g(x) x 4- f(x) x 2 + X(x,x) , 


'2 

where X (x,x) = o(x ) as x + 0, Equation (8) now takes the form 


x + f(x) x 2 + g(x) x + h(x) + X(xjX) = 0 . 


By analogy with the analysis of equation (4) we adopt the following form 
as a candidate for a Liapunov function, 

.2 

(9) V = — ex P [2F(x)l + H(x) , 


where 



f(s) ds 


g(x) + x 


x 

; H(x) = J h(s) exp [2F(s>] ds. 
o 


X(x,x)^| exp [2F(x)] 


The function (9) will be a Liapunov function for (8) if h(x) and 
g(x) are positive in some neighborhood of the origin; the extent of the 
region of asymptotic stability will depend on the nature of the 
function X(x,x). 
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Example. Rayleigh’s equation. 


*• • o # 

x - qp x +hx+x=0. 

Here h(x) = x, g(x) = p and f(x) = 0; thus 



* *2 / ' * 2 \ 

■V = - yX |1 - x j 

In this case V is the usual quadratic form. 

5. Third order equations. 

In this section we consider the third order equation 

(10) x + f(x) x* + g(x) + h(x) =0. 

Again, our aim is to put the new equation into 

* • 

treated equation. Define F(x), G(x), and H(x) 
equation (10) as 

(11) dt + + §( x ) + Mx) = 0 

The first two terms of (11) suggest the function 

V 1 = 2 D* + F( ^ 2 + G( * } * 

This, however takes no account of the third term h(x). We therefore, seek 


a form resembling a previously 
in the usual way, and write 



a Liapunov function in the form 


where is as above, and is to be determined so as to simplify the 

• • • • 

resulting expression for V. Since V = it is simpler to choose 

• • • 

V 2 so as to simplify V and further require that be given as a derivative. 

Differentiating V with respect to time along a trajectory of (10) 
yields 


V -['* + F(x)][ - g(x) - h(x)J + g(x) x + V, 
= - g(x)F(x) + |^V 2 - h(x) x - h(x) F(x)J 


• * 


The first term above, - g(x) F(x), corresponds to previous results and 
should be retained. 

• • 

In order to choose we assume that F(x) has a continuous first 
derivative; write 


F (x) = a x + . . . , 


and 

• ® 00 

F(x) h(x) = a x h(x) + [F(x) - ax] h(x) 

= a ~ [a H(x) ] + (f(x) - a x)h(x) . 

This suggests that we take 

V 2 = dT C x + a H ( x )1j 
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Then 


- h(x) x - h(x) F(x) = x 2 l/(x) + j^ax - F(x)J h(x) 

We thus, adapt as a candidate for a Liapunov functipn for equation (1) 

(12) V = y^x + F (x)J 2 + G(x) + x h(x) + a H(x) 

V = - [g(x) + h(x)J F(x) + x 2 h (x) + a x h(x) 

Example, Barabashin's equation. 

For the equation 

• •a •* • 

(13) x + a x + g(x) + h(x) = 0 , 

The above analysis gives the Liapunov function 

• 9 • 

V = - a x g(x) 

The V given here is the same as that obtained by Walker {2;p 87} using his 
modification of Infante’s integral method; it yields the following set of 
sufficient conditions for global asymptotic stability 

• .* 

.a g(x) - h (x) > 0 ; x ={= 0 

x 

a > 0 ; x h(x) > 0 ; x jO . 

' ’ - 16 - - • • 


x] 2 + G(x) + a H(x) + x h(x) 

+ x 2 h (x) . 



V -* 03 for x -»■ 00 


X X 

Ri'+.i ] 2 + [ g(s) ds + a ^ h(s) ds > x h(x) 


Schultz {5;p 39} and Puri {5;p 46} treat Barabashin's equation by 

the variable gradient method. Puri obtains the function 

* * 

x x 

V = ~ (x + a x)^ + a ^ h(s) ds + j g(s) ds + h (x) x 
P o o 

v - - * 2 [a *£> - h'(x)] , 

r X 

which yields almost the same conditions for global asymptotic stability as 
does the function V. 

6. Some other forms. '• 

In this section we present three examples which are not in any of the 
forms specified in the previous sections. These examples are intended to 
illustrate the gaps that develop: .as a method for generating Liapunov 

functions is extended. 

Example A. Walker { 2 ; p 89} studied the equation 

• • • • • • ’ 

(14) x+bx+(x + cx) =0, 


This equation is not of the form of equation (10), but bears some similarity 
to it. 
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Write (14) as 


* ® « 


— (x + bx) + (x + cx) m = 0 , 


and let 


V = i [x* + bx] 2 + feLJ— + v 

2 lX DXJ + C(m + 1) 2 * 


where 


(x 4 ex) 


m+1 


C(m + 1) 


= G(x) 


I 


(x + cx) dx 


Differentiating we get 


V = - (b - x (x + cx) m + V 2 


This form of V does not suggest an appropriate form for V^. If, however. 


we do not replace x by - bx - (x+cx) m in the evaluation of V we obtain 
the expression 


(15) 


• -I « • 1 • -I * * • • 

V s= (b - — ) x x + b (b ) x x + 

c c 2 


The order of x is now reduced by extracting an exact differential from 
the first two terms of (15). Thus, 

v.| t [; ;• + b ; z ] (b - i) ; 2 + v 2 


The choice of V„ is now clear. 


h ■ it r* * +bx2 3 - 1) 
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With this choice of we get 

v-K* +7 X ] 2 + h 2 (bc - 15 x2 + C 15T) < x + cx > m+1 • 
v--( b -|)i 2 • 

The Liapunov function obtained here is the same as that obtained by 
Walker. 

Ingwerson {7;p 39} and Szego {8 ;p 52} consider equation (14) for 
the case m = 3, and obtain very similar results. Ingwerson uses a linear 
analogy in his work while Szego uses a modification of Zubov's technique. 
Example B. Ingwerson {7;p 45} considers the third order system 

• t « • • • *2 

(16) x + b ± x + (b 2 + c 2 b 3 ) x + b 4 (x + c^x) + b 3 x = 0 . 

Let 

y = x + c 2 x 

With this notation equation (16) becomes 

IfL x + b x x] + b 2 X + b 3 y + b 4 y 2 = 0 . 

• • 

Since the coefficient of x in the above equation is unity, it would be more 

■ * » 

convenient to make the coefficient of x in y equal to unity. Thus, let 
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and our equation now becomes 


d_ 

dt 


C x + vO 


+ b 2 * + b 3 C 2 z + b 4 c 2 2 Z 2 


= 0 . 


We now take V of the form 


V = |-[x + b ± x] + b 2 


2 „ ,3 

'2 2 ' "4 "2 3 ~ ' '2 


7 o 7 

c„ — + b, C „ 2 — +v n . 


Differentiation yields 


(17) V - - b 1 b 2 x 2 - 




* r» • 

, . x , , 2 , x 

be x + + be x + — 

3 2 c 2 4 2 c 2 


] * + i 


We now look for an appropriate derivative in the second term of (17) 
to choose as V 2 » After some rearrangement of terms we find that 


[‘.-yiv (* 


k-},] t[ 


, b 4 C 2 • , 2b 4 C 2 , . , 

+ — : x +, x I I + 


b 3 C 2 


b 3 2 , b 4 3 
— x + — x 


■)J 


! ] • 


Choose 




b. - - 
1 c. 


d r b 3 2 . b 4 3 

L~ x + ~ x 


] : 


- 20 - 



Then 


v . - bl b 2 ; 2 - (b^-D. [b 3 ; 2 (i + 


V 2 • ^ 2b 4 

b 3 b 3 


)] 


•t • • *o ^0 *0 b«c~ ^ r 2 b / c T 

v.-i(x + v> + r x + -p z [i + -^ z 


(h-cjr x2 ( 1 + 3b7 x ) 


The above function gives the following sufficient conditions for 
asymptotic stability: 


b ± > 0 ; i = 1, 2, 3, 4 


b l c 2 =• 1 


1+ l£i (; + as)> 

b 3 V c 2 / 


1 + 


2b 4 C 2 


z > 0 


2b 4 

1 + — x >0 
3b 3 



Ingwerson claims global stability, but his Liapunov function is not 
given in the report , 
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Example C. Our last example is a fourth order equation considered by- 
Walker { 2 ; p 96} : *- 

8 * * « • • • • • « Q 

(18) x + 4x + 5x + 2x + cx = 0 . 


Even though equation (18) is of a relatively simple form, containing only 
one nonlinearity, a good deal of ingenuity is required to produce a Liapunov 
function. . * 

The terms in (18) may be grouped, and the equation rewritten as 


(x + x) + 3(x + x ) + 2 ( x +x)+ ex' 3 = 0 


This suggests that we let 


y = x + x 


and write the equation in the form 


(19) ^r( y + 3y) + 2y + cx 3 


= 0 


Equation (19) suggests that we try to construct a Liapunov function 
in the form , 


1 * * *2 “2 
V = “ ( y + 3y) + y + v 2 * 
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Differentiating, 


« 9 ry Q • • • •• • * • • *" , • 

V = - 6y - cx ( x + x + 3x + 3x ) + V 2 

2 ••• j r o »* • o * 

+ 3cx ( x x + 4x ) - c — j x (x + 4x) + — x J + . 


We thus take 

• d r 3 • , 3 41 

v 2 ■ c 3 ? L x x + 4x + 4 x J • 

which yields 



( 20 ) 


V=-6(x + x) + 3cx ( x 


*2 

4x Z ) 


The above form of V does not lead to satisfactory stability criteria; 
therefore we shall redefine , Equation (20) may be written as 


V = - 6(1-2 x ) 


(-S’- 


2 "2 21'* ** 189 "2 

12 cx x - — x x - —rr x 

64 Z 


Redefine as 


dt 


v, . c 4r[x 3 ( ;- + 4; )+ fx 4 ] + f ^ ( ; 2 > 


4 dt 


The last term of the new is designed to simplify the last three terms 

• • 

in the expression for V . With this definition of we get the final result 







II. An attempt was made by Dr. J, P. Clay to find out if any theory could 
be developed regarding the underlying structure of the class of all Liapunov 
functions (yielding asymptotic stability) admitted by a given system of 
differential equations. Coincident with this, the class of all systems of 
differential equations (of a given type) which admitted a given Liapunov 
function was studied. 

This attempt yielded only a few results of an elementary nature, and 
was abandoned as being premature. The generality of the approach was far 
too great to yield concrete results until more is known about specific 
classes of equations and their associated Liapunov functions. 

The following is a brief description of Dr. Clay's ideas. 

Let W be a closed neighborhood of the origin in n-dimensional 
Euclidean space, and let V(x) be a positive definite mapping of W into 
the reals having continuous first partial derivatives. 

We define 

Ly(W) = {f/f_: W -*■ E n ; _f (0) = 0; VV • f is negative definite} . 

Basically L^(W) represents the set of all autonomous systems of 
differential equations for which the origin is an asymptotically stable 
isolated singular point, and for which V is a Liapunov function. 

The following results hold: 

(1) f, g e L v (w) f + g e L V (W) 

(2) f s Ly(w) and a > 0 - — af e L V (W) 
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(3) Ly(w) CZ C(W; E n ) = the set of continuous functions 

on W into E n , 

(4) L y (w) is a convex topological semi-group of C(W; E n ). 

Let f_:W ->- E with f_(0) = 0 and f continuous. 

We define 

(W) = {V/V:W R; V has continuous first partial derivatives; 

V is positive definite; VV*f_ is negative definite}. 

• 

Assuming that the system x. = £(,2.) has an asymptotically stable 
singular point at the origin, then (W) is the set of all Liapunov 

m 

functions defined on W for the system x = f(x). 

The following results hold: 

(1) For any V^, V 2 e L f (W), and a > 0 

(i) V x + V 2 s L f (W) 

(ii) aV^ s L f (w) 

(iii) V x V 2 e L f (W) . 

(2) If V £ L^(W), then S £ (V) f*) L^ (W) =}= <{>. for any e > 0. Here 

S £ (V) = {V* / sup |V(x) - V*(x)[ < e } 
x 
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(3) (W) is a convex topological semi-ring of the ring of continuous 

functions on W into R, . 

An attempt was made to investigate the structure of L^(W) and L^(W) 
using topology, algebra, and functional analysis but no substantive results 
were obtained. ' 

The main body of this part of the report consists of our efforts to 
uncover some structure of a more specific class of differential equations. 
Specifically, we wished to classify in a meaningful way those systems of 
autonomous differential equations which admit of a quadratic form as a 
Liapunov function for the determination of the stability or instability of 
the origin. 

We considered the system 

• * 

(a) x - f(x); f(0) = 0 , 

where f.(x) is an analytic function of x. 

It is well known that if (a) has the form 

(b) x = Ax + F(x) , 

where A is a constant nxn matrix and is analytic with no terms 
of degree less than two, then (a) admits of a quadratic form as a Liapunov 
function if the eigenvalues of A either all have negative real parts or at 
■ least one eigenvalue has a positive real part. We conjectured that if the 
matrix A has no eigenvalue with a positive real part, and also no eigenvalue 
of multiplicity greater than one with zero real part, then (b) admits of a 
quadratic Liapunov function. This conjecture turned out to be false as the 
following example shows. 



x = y 

} ' + 4 

y = - y + x . 

For the system (c) 



with eigenvalues A^ = 0, = -1. However, for' any quadratic form 

V =| (ax“ + 2bxy + cy^) 

we get 

* 2 S 4 

V = (b-c) y + (a-b) xy + bx + cyx 

• 5 

Along the line y = 0, V = bx which may take on both positive and negative 
values in any neighborhood of the origin. 

On the other hand, one can produce any number , of examples of equations 
which admit of a quadratic Liapunov function (or, in fact, of any given 
Liapunov function) by the following device. 

Let V(x) be a given positive definite function, and let P(x) be any 
sign constant scalar function. Then if G(x) is such that C? (x) * W(x) = 0 
the equation 

e * 

(d) x = P(x) W + G.(x) 
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has V ( x) as a Liapunov function. The characterization given in (d) is, of 
course, not in a useful form except to produce examples. 

The next section of this report describes our efforts to broaden the 
scope of the problem in order to achieve meaningful results. 
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Equations with homogeneous right hand sides 

1. Introduction, 

We devoted considerable effort to systems of the form 

(1) x ± = b ± (x); i = 1, . . . ,n, 

where . . . ,b^ are homogeneous polynomials as a first step toward 

finding conditions on the more general system 

(2) . x ± = f±(x); i = 1, . . . ,n 

which imply the existence of a homogeneous polynomial Liapunov function. 
Suppose the right side of (2) is analytic in ). Then 

CO 

= I f ^ . (x) » i = l»...jn, 

j=l J 

where f (x) is a homogeneous polynomial of degree j . 

If fil*0 for some i, then (2) has a linear part; this is the classical 
situation studied by Liapunov. Suppose k > 1 is the smallest integer for 
which f.^0 for some i. If k is even, (2) is unstable. [I;p96], Hence we 
assume that k = 2m + 1, in which case (2) can be written 

(3) x ± = b i^— ^ + s i^-) ! 1 ~ I» • • • 

where b_^ is a homogeneous polynomial of degree 2m + 1, and g. contains only 
higher degree terms. If V (x) is a Liapunov function for (1) (V positive 
definite, and V < 0) it is also a Liapunov function for (3); hence our 
interest in (1) . 

It follows immediately from homogeneity [l;p97] that if (1) is 
asymptotically stable, it is asymptotically stable in the large; thus we 
shall always refer merely to asymptotic stability. 



As a starting point for the investigation of the existence of homo- 
geneous Liapunov functions for (1) we considered generalizing the following 
theorem of Liapunov [2] : 

n 

Theorem 1, If the eigenvalues of the matrix A = (a..). . _ 

' ij = i 

have negative real parts and if U is a negative definite homogeneous form 
of degree 2k, then there exists a positive definite homogeneous form Q of 
degree 2k such that ’ 


I Cl a . x ) = U. 

i=l i=l J 3x. 

J r 

This theorem guarantees the existence of Liapunov functions (and 
asymptotic stability) for the linear system 


n 

(5) x. = J a., x. ; i=l,...,n. 

i y 3 

An analogue for the nonlinear system (1) would give conditions on the 
functions b^,.,.,b n such that, for a given negative definite form U of 
degree 2k + 2m, there is a positive definite homogeneous form Q of degree 
2k that satisfies - 


(6) ' I b (x) =U(x) 

i=l 8x. 

x 


This question, however, is not well posed unless we restrict 13 ; if V is 
of degree 2m + 2k it has 

2m + 2k + n - \ 
n - 1 
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coefficients, while, if Q is of degree 2k it has only 
/2k + n - l\ 

A n - 1 / 

coefficients. Hence for m > 0, the dimension of the space of homogeneous 
polynomials which can be written in the form of the left side of (6) is 
lower than the dimension of the space of homogeneous U's of degree 
2m + 2k. Thus it is necessary to add the explicit assumption that U lies 
in the lower dimensional subspace of homogeneous polynomials which can 
be written in the form (6) for any homogeneous Q of degree 2k. 

2. The eigen-polynomial problem. 

The hypothesis of theorem 1 on the eigenvalues of A is necessary and 
sufficient for the origin to be asymptotically stable for (5). The 
classical notion of eigenvalue can be generalized for (1) by calling a 
homogeneous polynomial A(x.) of degree 2m an eigen-polynomial of the set 

} of index 2k if there exists a non-trivial homogeneous polynomial 
P (x) of degree 2k + 1 such that 

n 

(7) l b.(x) 8P = A(x)P(x). 

i=l 1 3x. 

l 

We call P an eigen-form of (b^,...,b }. (It is not necessary for 
{b^,...,b n } and P to be odd and A even to make good sense algebraically, 
but we restrict our attention to this case because it is the only one of 
interest in stability theory.) 

This is a generalization of the classical eigenvalue problem: if 

k = m = 0, A reduces to a constant, and • 



u 

b.(x) = 7 b.. x. 

I s - An ij 3 


3=1 
n . 


p (x) = l a A X -i 
3=1 3 3 ’ 


where b.. and a. are constants. Now (7) reduces to 
13 3 


n n n 

Y a. f 7 b,. x.i = A I a. x. 

ML ^1=1 13 V 4l 3 3 * 


n Iji 

Thus X is an eigenvalue of B = (b„) and A = (a^,...a n ) an 

13 i j =1 

T 

eigenvector of B . 

The relevance of the above discussion to systems of the form (1) is 
as follows: 

If x (t) is a trajectory of (1), relation (7) implies that 


d P(x(t)) = X(x(t)) P(x(t)); 

dt 


hence 


v- 

P(x(t)) = P(x(0)) expj^J* X(x(t) dsj 


So that if X is negative definite, P(x) approaches zero along any trajectory 
of (1). 

As far as we know there is no well developed theory of eigen-polynomials 
and eigen-forms as defined here. We have confined our attention to the case 
k = 0, for which (7) takes the form 


( 8 ) 


n r n 

l a.b (x) .= X (x){ l cc x ) 
i=l 1 1 v i=l 1 x / » 


- 33 - 



with constants a^,...,a . The eigen-form 


* -l “i x i 

i=l 


satisfies the equation 

y = X (x (t ) ) y 
along a trajectory of (1). 

Suppose now that there are n linearly independent eigen-forms 


n 

y. = T a., x. ; i = 1, . . . ,n 
i ij 3 


(which implies that the matrix A = (a„) is non-singular) with associated 
eigen-polynomials A^, A^j***?^* then (1) is equivalent to 


* _1 

y i = X i y i ; 1 = 


where ^(t) - (y^ (t ) , . . . y^ (t ) ) , 


This is the analogue of the diagonal form for (5) . The question of the 
existence of n linearly independent eigen-forms, or more appropriately, the 
determination of all canonical forms for (1) , requires an analogue of the 
Jordan canonical form for matrices. We have not attempted to develop such 
a theory because we thought it wise to see first if we could settle the 
question of stability for the diagonal system. Our efforts in this direction 
are reported in sections 4 through 8. However, we did obtain an algebraic 
formulation of the eigen-polynomial problem for ,k=0, which we report here 
for completeness. 
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a) 


X 1 ^ x l J °) < 0; X 2 < 0 

b) at every (x^jX^) 4 s (0,0) at least one of X^ (x^jX^) 

and ^2 (x^,x 2 ) is negative. 


Proof : The proof makes use of the basic relation (15) which holds along 

any trajectory of (14). 

05) _i. (X, - X 2 ) - (A 1 - A,) x, +- Jh xj 

This relation is obtained by applying Euler’s identity for homogeneous 
functions : 


3X x + 3X x = 2mX 
a X;L 3 x 2 

to the derivative of X^ - X^ along a trajectory. It is clear from (15) 
that a line along which X^ = X^ is a trajectory of (14). Moreover (15) 
shows that X^ - X^ is equal to its initial value times an exponential, 
and hence X^ - X 2 does not change sign along any trajectory. 

We prove first that every trajectory of (14) approaches the origin 
(quasi-asymptotic stability). It follows immediately from the representa- 
tion (13) and hypothesis (a) that any trajectory which starts on one of the 
coordinate axes must approach the origin, and in fact, monotonically . Now 
assume that x^(0) 4 s 0 and ^(O) 'f 1 0. Since X^ - X^ does not change sign 
along a trajectory we may assume without loss of generality that 


x 2 -x 1 >0. 


From (13) 


x 2 (t) _ X 2 (0) 


x x (t) 


x 1 (0) 


“ P [J 


{X 2 (x (s)) 


- X 1 (x (s)) } ds^j 
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Def ine 


so that — approaches a limit l (possibly infinite, but not zero). 
X 1 



r r 


Then and ^ approach limits : 

lim E, = 1 = /T hT" . 
t “ . • 

<16) . : ' 

lim = C 9 = / 1+2. 2 . 

t . ' 

Now, (^-^9^2) ~ 0 a n d it follows from hypothesis (b) that 

If ~ a w ^ ere a > 0* Thus there exists a T such that t >_ T implies. 

&7) h ( *L’V < - f- r 2m 

From (13) 

t 

[x 1 (t)| = J X;L (T)| exp LI A, (x(s)) ds"| 

■ T 

-t „ 

<_ |x 1 (T)| exp 




t 2m ~ 

<_ (T) 1 exp | 

r o' 

_ " 2 

f ( X;L (s)) ds J 


T 


From (17), jx^j is decreasing for t _> T; therefore it has a limit | x^j . 
The assumption that | x^J=}=0, along with the last inequality of (18), leads 

to a contradiction. Hence lim x. (t) = 0. Since C-,4 5 0 (because l =j=' 0 in 

X 1 t - 9 - » 1 

(16)- and r = — , lim r (t) = 0, which completes the first part of the 

4 1 t ->0 


[- f- i 


zm 

(r(s)) ds 


] 


proof. 



The above proof shows that once a trajectory of (14) enters a portion 
of the plane in which both X^ and X ^ are negative it remains there and 
r(t) "*■ 0 monotonically, Thus, starting at a point where r(x^(0), X2(0)) = r Q 

and say X^(x^(0) jX^CO)) < 0 the maximum value r can assume is less than 
2 2 

*2 (0) + n where n is the number of smallest absolute value such that 

2 2 

^2 - 0. Hypothesis (a) now assures us that x^ (0) + n can be 

made as small as desired by choosing r Q sufficiently small. This completes 
the proof of the theorem. 


5. A conjecture on the diagonal case for general n. 

We sought to generalize theorem 2 with the following conjecture: 
Conjecture 1. The system 

x ± = X ± (x) x ± ; i = 1, 2,.. .,n. 


where X^j...^^ are real homogeneous polynomials of degree 2m, is 

asymptotically stable if and only if for any partition [i^, . . . , i^J » 

[i, i ] of N = [l,...,n] into disjoint subsets, (either of which may 
xC* x n 

be empty) at least one of the polynomials X. ,...,A, is negative at any 

1 1 H 


non-zero point on the hyperplane x. 


= x. 


• « « *”■ X. 

1 


o. 


k+1 A k+2 n 

We expended considerable effort to prove this conjecture, but without 
success. At this point we can only surmise that either the hypothesis is 
too weak, or the theorem is very deep. The crucial defect in our attempts 
to prove it seems to be the lack of a generalization of the relation (15) 
which implied that X^ - X 2 is sign constant along a trajectory. 

We have obtained some partial results which are presented here. 

In the following discussion we may assume that no trajectory intersects 



a coordinate hyper-plane, for if x (t ) = 0 for some t >0, then 

no o — • 

x n (t) = 0 for t >_ t Q , and (19) reduces to a system of n - 1 equations 

which satisfy the hyperthesis of the conjecture with n replaced by n - 1. 

Lemma 1 , Let the system (19) satisfy the hypothesis of conjecture 1, 

and let x(t) be a trajectory along which one of the A^ dominates the others 

for sufficiently large t. Then lim x(t) =0. 

t-x» 

Proof : Without loss of generality, assume that 

A n (x(t)) - A^ (x(t) ) >_ 0 ; t >_ T, i = l,...,n-l. 

From (13), if t > T 


x (t) x (T) 
n _ n exp 

X i (t) ~ x i^ 


£ j U n (2£(s)) ~ A. (x(s))}J 


ds 


hence 


n 


, being non-decreasing for t >_ T, approaches a limit l. 

”X r " X 


(possibly infinite, but not zero). Let 


x - (x 2 + . . . + x 2 ) 
l n 


1/2 


and 


Then 


^ — x , i — 1, ... . ,n. 


lim £. 


. (t) = 5 = _1 _ ( 1 +1 + ... -fr 1 +1} 

1 1 1 Vh 2 V / 


- 1/2 


where i =1. Suppose 4 = 0 for i in the set [i, .....i, ] 
n x 1 1 k 

2 

i in [i^ + ^, . . . ,i^] . The first set is not empty (since 
and by hypothesis contains an integer j such that Aj(£^,.. 
As in (18) , there is a > T such that for t >_ T^ 


and = 0 for 
x 

+ ... + 5 2 = 1 ) 
n 

,5 n ) = -o (a >0) 


9 



( 20 ) 


2 


x. (t)l <_ } x. (T )j exp 
3 3 ^ 


i 


y» ’i 


Lim | x . ( t ) 1 exists because |x.(t)| is decreasing for t >_ T^, and the 

£-Ko 0 d 

assumption that this limit is not zero leads to a contradiction via (20) . 


Thus lim x.(t) = 0, and since r = 
t-*» ^ 


_ 1 




and f 0, it follows that 


lim r(t) = 0. This completes the proof of the lemma. 

t-*» 


Corollary 1. If X^,...,X n satisfy the hypothesis of conjecture 1 
and X^ - X^ is positive semi-definite for 1 <_ i <_ n-1, then the system (19) 
is asymptotically stable. 

Corollary 2. If X^,..., X^ satisfy the hypothesis of conjecture 1 
and if x(t) is a trajectory for which approach limits, then 

lim x(t) =0, 


If lim % ~ £ exists for i = l,...,n, Corollary 2 implies that 

^1* * • * 9 — 0 ; 1 • • • jH* 

Thus the hypothesis of conjecture 1 and the assumption that £^,..., E, n 
have limits are incompatable, unless the former implies the existence of 
at least one point where X^,...,X are all non-positive. This is true for 
n = 2, but we have not been able to show it for n > 2. In conjecture 2 
(Section 8) we strengthen our hypothesis so as to remove this difficulty. 
However, the stronger hypothesis brought us no closer to a solution. 

6. An estimate. 

In trying to translate the hypothesis of conjecture 1 into usable form, 
we discovered the follox^ing lemma. 
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Lemma 2. Let A^,...,A n satisfy the hypothesis of conjecture 1; let 
2 , % v 2 


S = and N-S = ti k+l , *** ,i n‘ 1; let 


■S <S> ' J *1 
leS 


and' 


r 2 (x) 
N-S 


1 


ieN-S 


For each S =f= N, there exists a constant A(S) such that if 

( 21 ) 


then 

(22) 


A^(x) >_ 0 ; ieS 


r 2 OS) Ji A ( s ) r 2 (x) . 

S N-S 


Proof. If (22) is false for some set S, then for each integer j 
there exists a. point x^ which satisfies (21) and for which 

2 ,,(j). 


• s (x ) > j r (x ) . 


2 (x (j> ' 

N-S 


From the 

homogeneity of A., x^ 

(23) 

r s 2 (x Cj) ) - 1 

Then 


(24) 

r 2 (x^) < 


N-S 

and 



1 

3 


r '<x vj O = r (x VJ ') + r 2 (x^) < 

N-S 


1+ J 5 


.(j) 


A„Ci>x _ _ 2 ,„(j), 

i, > 

hence {x vjy } is a hounded sequence and has a limit point x . From (23) 

r c (x) = an< i from (24), r (x) = 0. Therefore x is a non-zero point 

N-S , 

on the hyperplane defined by x = 0 for i e N-S. From the hypothesis of 



conjecture 1, A^x) < 0 for at least one i £ S. This is a contradiction; 

A 

the definition of x. and the continuity of imply that A_^ (x) >_ 0 . This 
proves lemma 2. 

Lemma 2 at first seemed to he quite promising. Taking S to be the set 
of all i for which A >_ 0 at any given point, and N-S the non-empty set of 
all i for which A^ < 0, it can be seen that (22) gives an upper bound for 
the sum of the squares of those coordinates whose absolute values are in- 
creasing, in terms of a similar sum of terms whose absolute values are 
decreasing. In fact (22) can be replaced by 

r 2 (x) < (1 + A(S) ) r 2 (x) . 

N-S 

Over any period of time during which no A^changes sign, this gives a 
decreasing upper bound for the distance from the trajectory to the origin. 

We have not been able to use this seemingly favorable state of affairs 
to prove the conjecture. 


7. Two theorems. 

In attempting an induction proof of conjecture 1 we were led to write 
(19) in the form 


(25) 


— (A^( x ^> • • • » x n -l,0) 1> * * • »n-l a 


x = A (x) x. , 
n n— 1 ’ 


where <Jk is homogeneous of degree 2m-l. 

Theorem 3. Let the system 
(26) y^ = A_^ (y^, • • • ^ ^ ~ l,,..,n— 1 , 

be asymptotically stable. Then there is a constant 0 < 6 < 1 such that if 



|g (t) I < 6 along a trajectory of (25), then lim x(t) =0. • 

n — . • — 

Proof, We use a. device due to Zubov. (1}. The asymptotic stability of 

(26) iaplies that any trajectory j(,t) of (26) satisfies 

(27) . |x(t) ] = 0(t" 1/2m ) as t-*» , : 

(Krasonskii [3]). • . ' 

* * / 

Zubov showed that this guarantees the existence of a Liapunov function 

V(y 1> ...,y n _ 1 ), continuously differentiable with respect to (y^,...,y n 
with derivative along a trajectory of (26) given by 


(28) 


dV 

dt 


n-1 

I 3V . , 

1-1 17 . x ± ^v-’ y n-v 0) 


. 2 , , 2 m+k 

- ( yi + • - + vd 


where k is a sufficiently large positive integer. 
Explicitly, • 

V ( 


u > = J j i (t.u) 1 2mrt2k dt , 


where (t, u) is the solution of (26) for which £ (o, u) = _u = (u^, , . . » u n _^) , 

and [ | is the Eucledian norm. The existence of the integral follows 

from (27). V can be shown to be positive homogeneous of order 2k, 

Suppose x^(t),...,x n j(t) are taken to be the first n-1 components of 
a trajectory x(t) of (25); let 

if>(t) = V(x 1 (t),...,x n> _ :L (t)), 

2 2 2 
p = x 1 +...+X 


and 


n-1 


2 2 2 2 2 
r - x. -K , ,+x = p + x 

1 n n 


From (25) and (28) , 
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* 2m+2k 

t - - P 


+ x 

n 



3V 

9x 


x . 4>_- * 

i i 


There is a constant M > 0 such that the second term on the right is 
dominated by 


hence 

* 2m+2k .. | _ i/T f 2 N m+kv 

1 1> < ~ p (1. - M | C n i (1 - C n ) ) . 

choose 6 so that I £ I < 6 implies that 

n — 

2 ' 

M | C n | (1 - 5 ) <_ 1 - a (0 < a < 1, arbitrary) ; 

p 

then if | £ | <_ S for all t, 

t 

(29) . .*(t) < 4(0) - a \ p(s) 2m+2k ds. 

■ ° . 

■ • j 

and is strictly decreasing. Let lim ^(t) = \p ; we xrant to show that 

t-Xo ■ ■ 

i[i = 0. If tj; =j= 0, then p is bounded away from zero, because for all t, 

* _• - 2k 

0 < i{/ <_ if»(t) = V(x 1 (t),...,x n _ 1 (t)) <_ p p (t), 

where p = min V(x^, . . . ,x n _^) on p = 1. Substituting in (29) yields 

(30) *Kt) < <K0) - 1+ £ . 

a contradiction, since it implies that ip( t) < 0 eventually. Since $ = 0, . 
p ss 0 because V is positive definite; therefore lim |x j = 0 also (if not 

t-no 

then lim | £ | = 1, violating the assumption that j £ | <_ 6 <1 for all t) . 

t-*» 

This completes the proof of theorem 3. 

The hypotheses of theorem 3 can be weakened to require that 


M 


n 


a - c ) 


2 N — (m+k) 


2m+2k 
p ; 



| | | _^_-<5 for all sufficiently large t. • 

Theorem 4. Let the system (26) be asymptotically stable and j£(t) 

be a trajectory of (25) for which lim x (t) = 0. Then lim x(t) = 0, 

t -* 00 ^ 

Proof, Using the notation of theorem 3, there exists a 6 in (0,1) 
such that either | £ j > 8 or 


(31) 


o 

ip < 


2m+2k 

a p 


For any t such that (31) does not hold, j ^ j > 6 and 

rxo\ ir+\ /' 2k 2k /'i X n \ ) 2k 

(32) f(t) < n p <_ nr _< ri — 1 J , 


where n = max V(x^, . . . ,x n _^) on the sphere p = 1. If either (31) or (32) 
holds for all sufficiently large t, we are finished; if this is not the 
case, there exists a sequence t 0 <t j_ <X 2 < *** a PPtoaching infinity such 
that (32) holds for < t < T 2j+1 ’ and holds for 


T 2j+1 - Z - T 2j+2 . Then 


4>(t) 


- n 2k , 


T 2j - t - T 2j+1 


and 


*(t) 1 ^^ T 2 j+i ) 1 11 




X 2j+1 - t - x 2j+2 


Since x^ approaches zero, so does ip , and therefore p approaches zero. 
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This completes the proof of theorem 4. • . 

Corollary. If the system (26) is asymptotically stable and X^ is 
negative definite, then the origin is asymptotically stable for the 
system (25). 

8, A stronger hypothesis. . . 

The hypothesis of the following conjecture is stronger than that of 

conjecture 1; if has the advantage of guaranteeing the existence of a 

region where are all negative. 

Conjecture 2. Let A^ s A 2 ,...,A n be a set of constants such that 

2 2 

A. > 0 and A- +...+A < 1. and let X. < 0 whenever jc.l > A.. Then 

x — I n .x x x 

the system 

• 

x. = X. (x) x. ; i = I,. . . ,n 
x x ■— x 

is asymptotically stable. 

We have not been able to obtain sharper results with this stronger 
hypothesis. 
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